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Abstract 



We investigate constraints for including bulk and brane matter in the Randall- 
Sundrum model. In static configurations with two zero thickness branes, we find 
that no realistic brane matter is possible. We also consider the possibility that the 
radion has stabilized by dissipating its energy into the bulk in the form of some 
unspecified matter, and find the Randall-Sundrum cosmological solutions in the 
presence of bulk ideal fluid. We discover that the metric is necessarily in a static 
configuration. We also discover that there is only one allowed equation of state for 
the bulk fluid, p = p, corresponding to the stiff ideal fluid. We find the corresponding 
brane cosmologies and compare them with the Friedmann- Robertson- Walker model. 
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1 Introduction 



Randall and Sundrum have recently suggested a novel solution to the hierarchy problem 
involving extra dimensions |1], 0, which has attracted much attention. The RS proposal 
assumes a five-dimensional spacetime where the extra dimension is a S 1 /Z 2 orbifold, 
with branes located at the two orbifold fixed points. The branes are thus at the spatial 
boundaries of the bulk spacetime. The brane tensions together with a negative bulk 
vacuum energy provide the source for the 5-d Einstein equation, which yields as the 
solution the metric 

3 

ds 2 = e" 2%l (-dt 2 + E^) 2 ) + dy 2 . (1) 

i=i 

The exponential warp factor rescales the physical masses on the second "TeV brane" , 
giving rise to a large suppression factor with respect to the Planck mass, and hence a 
possible solution to the hierarchy problem. However, this is achieved at the expense of 
fine tuning: the cosmological constants on the two branes must be of equal magnitude 
but with opposite signs. Also, the size of the extra dimension, or the radion, is a four 
dimensional modulus field whose value needs to be fixed at the right scale in order for the 
RS scenario to work. 

The RS model has been refined in many ways since its inception, but the question of 
radion stabilization remains unsolved. According to commonly accepted wisdom, moving 
branes would give rise to particle masses which change in time [[|, so that the radion 
should be very nearly stabilized by the beginning of nucleosynthesis so as not to conflict 
with observation. One possibility was proposed by Goldberger and Wise, who assumed 
the existence of a massive bulk scalar field with self- interactions on the branes M; see also 
||. [|. By integrating over the fifth dimension one then generates an effective potential 
for the radion, which has a non-trivial minimum. As the radion field settles into the 
minimum, the size of the extra dimension gets fixed. 

Using a bulk scalar field to stabilize the radion has also been studied in 0. Other ideas 
for fixing the size of the fifth dimension include gaugino condensation in a supersymmetric 
setting || and the possibility that the Hubble redshift might damp the radion so that it 
is almost stabilized ||f]. A phenomenological mechanism for stabilizing the radion was 
considered in ||. 

The implications of radion stabilization independent of the specifics of the stabilization 



mechanism were considered in [11, 12]. It was found that a fixed size for the fifth dimension 



requires a certain non-constant form for the yy- component of the bulk energy-momentum 



4 In fl(i|| the interesting possibility that an asymptotically constant radion emerges from the Einstein 
equation naturally without the introduction of extra degrees of freedom was considered. Unfortunately 
the Ansatz used actually precludes any dynamics for the size of the fifth dimension. 
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tensor. This form turned out to be determined by a constraint that related the yy- 
component of the stress tensor to its trace. This constraint could be understood in terms 
of the backreaction of the radion to the inclusion of matter. 

Whatever the actual mechanism for stabilizing the radion, it is natural to assume that 
the fixed size of the extra dimension results from dynamical evolution. In the Goldberger- 
Wise scenario, the radion would evolve from its initial state towards the minimum by 
dissipating the extra energy away. The possibility of the radion stabilizing via decay into 
Standard Model particles on the TeV brane was considered in j3|. However, it does not 
seem farfetched to assume that there are other fields in the bulk, which can couple to the 
radion, whether directly, via the GW field or via some other mechanism, providing an 
alternative channel for energy dissipation. Such additional fields could e.g. be a part of 
the field content of a supergravity theory in the bulk. The asymptotic final state in the 
bulk would then involve the radion field with a constant value dictated by the minimum 
of the effective potential, together with some bulk matter. 

If bulk matter exists, at large times one would expect it to be in the state of maximum 
entropy. This means that e.g. viscous flow should eventually get damped away, and that 
the bulk matter can be considered an ideal fluid. In this case, stress tensor of the bulk fluid 
will be spatially isotropic, in particular the ^-component will not differ from the other 
diagonal spatial components, unlike in previous investigations |7|, [10], [II], [12], |I|, [14], [15|, [16 



In section 2 we find the general form of the metric assuming homogeneity and isotropy 
with respect to the three visible spatial dimensions, one static extra dimension and no 
matter flow along this extra dimension. We then focus on two-brane systems with zero 
thickness branes, and make a distinction between static and non-static configurations. 
We investigate constraints on brane matter in two-brane systems, and show that in static 
configurations the only allowed forms of brane matter are a cosmological constant and 
domain walls moving at the speed of light. Single-bvane configurations easily allow generic 
matter, but we are interested in the possibility of adding another brane necessary for 
generating the hierarchy. We then comment on the scenarios for including realistic brane 
matter in non-static two-brane configurations. In section 3 we consider the case of an ideal 
fluid plus a cosmological constant in the bulk and find the exact cosmological solutions. 
Curiously, we find that the metric must be in a static configuration, leaving no room for 
realistic brane matter, and that there is only one allowed equation of state for the bulk 
fluid, p = p, corresponding to the stiff ideal fluid. We note that it is possible for both 
branes to have a positive cosmological constant, and one brane may even have a zero 
cosmological constant, a result discovered in |17[] and emphasized in 0. In section 4 we 
discuss our results. 
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2 Warped spacetimes bounded by two branes 



2.1 General form of the metric 



The object of our interest is the Einstein equation in 4+1 dimensions, 



Gab = k Tab ■ 



(2) 



In the above, k 2 = 1/M 3 , where M is the Planck scale in five dimensions. The indices 
A and B run through time t, three spatial coordinates x l with infinite range and a com- 
pact spatial coordinate y. We use the standard notational convention and denote partial 
derivatives with respect to t and y by dots and primes, respectively. 

We begin with two general assumptions regarding the metric and the matter content: 

1. The spacetime is homogeneous and isotropic with respect to the three spatial coor- 
dinates x\ 

2. There is no matter flow in the y-direction. 

Assumption 1. The most general metric obeying homogeneity and isotropy with re- 
spect to the spatial coordinates x % is 



ds 2 = -nit, yfdt 2 + 2 c(t, y) dt dy + b{t, yfdy 2 + K,yj ^(dx*) 2 , (3) 



where r 2 = Y^=i( xl ) 2 an d 6K is the constant three-dimensional spatial curvature. While 
the general form of the metric includes a dtdy cross term, most of the literature has 
adopted Gaussian normal coordinates, where the cross term is absent (c = 0). The same 
assumption will be adopted in the present paper. We call the metric with c = the 
cosmological Randall-Sundrum metric in contrast to the original static Randall-Sundrum 
solution. For simplicity of notation, we put K = for the rest of this work. This does 
not affect our essential results. 

The Einstein Equation. The nontrivial components of the Einstein equation for the 
metric @ in Gaussian normal coordinates, with K = 0, are 



3 
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(4) 



We are interested in models with two parallel branes at the endpoints of the y- 
coordinate. In particular, we demand that the proper transverse distance between the 
branes is independent of time. So we make one more assumption: 

3. In the coordinates where the branes are at y = yi, y 2 , we have 6 = 0. 

Assumptions 2 and 3. No matter flow in the ^/-direction means that T ty = 0. Inserting 
this into the relevant component of the Einstein equation (H), along with b = 0, we obtain 
the result 

a(t,y) = A(t)n(t,y) , (5) 

where A is an unknown function. We are interested in cosmological solutions, and so will 
throughout the paper assume that a ^ 0. 

2.2 Motivation for a static fifth dimension 

The assumption of a static fifth dimension appears in most papers on the RS model. 
Since it does not seem obvious that a non-static fifth dimension would in general lead to 
conflict with observation, we present a brief discussion to motivate the assumption 6 = 0. 

The RS model was first envisaged as a solution to the hierarchy problem: the (four- 
dimensional) length scale increases as one moves from one brane to the other, which 
induces a change in the mass scales between the branes. The hierarchy between the 
Planck scale on the (hidden) brane at y — y\ and the TeV scale on the (observable) brane 
a t V — 2/2 requires that the functions n(t,y) and a(t,y) satisfy 



n(t, yi) 
n(t,y 2 ) 
a(t,yi) 
a(t,y 2 ) 



N 

Nh(t) 



(6) 



where N ~ 10 16 is the ratio of the Planck and TeV scales and h(t) is some function of 
time (to allow for the possibility of different cosmological expansion factors on the two 
branes). This is a rather weak constraint, since it makes no reference to behaviour away 
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from the branes. Therefore, in the literature there often appears a stronger and hence 
more pragmatic constraint. One makes the additional Ansatz that the condition (|6]) holds 
not only for the particular value y = y±, but for all values of y, with the value N replaced 
by some function f(y). In other words, the Ansatz says that the Einstein equation implies 
no preferred brane positions. This Ansatz implies that 

n(t,y) = n (t)f(y) 
a(t,y) = a (t)f(y) 

b(t,y) = b(t,y) . (7) 

Substituting (|7|) into the ty-component of the Einstein equation (f|) and assuming 
T ty = 0, we obtain the condition 

fh = . (8) 

For a nontrivial warp factor, the above equation can only be satisfied if b = 0. One 
can then redefine the coordinates t and y to set the functions n (t) and b(y) to unity, so 
that the metric reads 

n(t,y) = f(y) 
a(t,y) = a (t)f(y) 

Kt,y) = 1. (9) 

As an aside, we note that the redefinition of the y-coordinate to obtain 6=1 con- 
tains a possible problem: if the original function b(y) has zeros, the required coordinate 
transformation may be singular. In particular, this singularity may map either or both of 
the boundaries of the ^/-coordinate from finite values to infinity, making the fifth dimen- 



sion non-compact. This issue was discussed in [pL8| , and singularities in the RS scenario 



have more generally been considered in [17], [19], |2(J. In the present work, we simply as- 



sume that the metric is non-singular so that the fifth dimension remains compact in such 
transfer mat ions . 

We see that the factorized metric (§) and the condition 6 = can be motivated by 
requiring Randall-Sundrum-type solutions. If there is a weak time-dependence in b, we 



might expect to approximately recover the f actor izable metric (Bj), as noted in [21|. It is 



known that the metric (H) does not allow for brane matter (other than the cosmological 



constant) [15, 21]. Next, we discuss the constraints for adding brane matter. 



2.3 Constraints on brane matter 

In addition to presenting a mechanism to generate the hierarchy, the RS model had 
another remarkable feature: the localization of gravity into the vicinity of the Planck 
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brane. Clearly, one would like to preserve both of these features while adding matter onto 
the TeV brane. Are there then any constraints on possible equations of state for brane 
matter, and what features of the model are responsible for them? We will next focus on 
isolating the potential sources for trouble and find out what features the model should 
have in order to accommodate generic brane matter. 

The localization of gravity is easiest to achieve in the orbifold geometry where the 
branes have a zero thickness^]. This assumption will be adopted in the present paper: 

4. Both branes have zero thickness. 



Assumption 4. According to the above assumption, the stress tensor of brane matter 
is proportional to a delta function. Thus the Einstein tensor must contain delta functions 
as well. The orbifolding of the y-coordinate introduces discontinuities in the first y- 
derivatives of the metric, and delta functions in the second y-derivatives. Since the metric 
is assumed to be continuous, the only possible delta function contributions to the Einstein 
tensor come from these second ^/-derivatives, a" and n" . This severely restricts the input 
of the metric to the brane stress tensor (or vice versa). 
The brane part of the Einstein equation (H) reads 



3± a - 
6 V a 

1 fa" n"\ 

V 2 - + - 

cr \ a n 



k 2 T* 



k 2 T\ 



brane 



brane 



(10) 



The notations 5 and brane refer to the delta function parts of the derivatives and the 
stress tensor, respectively. With the assumption of spatial homogeneity and isotropy the 
general form of the brane stress tensor is 



B 



brane 



E Ul 7. ™' <&&g(-pm if) , p m (t) , p m (t) ,p m (t),0) 



m=l,2 



where the index m enumerates the branes, and we take V2 > Vi- On the other hand, the 
delta function part of the metric is related to the jumps of the first derivatives of the 



metric, which can in turn be expressed in terms of the continuous part of the metric [|22 

a "\s = 8(y-y m )[a'] 

m=l,2 

= E S(y-y m )(-ir +1 2a' c , (12) 



m=l,2 



5 Branes with finite thickness have been considered in |Tl], 12, [l9| 
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where [a'} is the discontinuity of a', [a'(y)] := lim £ ^ (a'(y + e) — a'(y — e)), and a' c is the 
continuous part of a'. An identical relation holds for n. Putting together (|10D, (pTT|) and 
fH2T), we have 



±2< 
6 a 

6 n 



y=y m 



y=ym 



-yPm (13) 
y(3p m + 2p m ). (14) 



The static configuration of branes. In the simplest scenario to generate the hierarchy 
between the branes, one assumed that the warp factor n is independent of the cosmic time 
t, so that n = f(y). In fact, there is always a coordinate system where this is true. Any 
two-dimensional Riemannian manifold is conformally flat, so we can change to conformal 
coordinates in the (t, |/)-subspace. In conformal coordinates, the metric (|3|) reads 

3 

ds 2 = -bit, yfdt 2 + b{t, yfdy 2 + a(t, yf • ( 15 ) 

We then orbifold in the ^-direction, so that we have two parallel branes at the fixed 
points at yi,y2- Introducing the assumption db/dt = 0, the metric reduces to 

3 

ds 2 = -~b{y) 2 dt 2 + biyfdy 2 + a(t, y) 2 ]T(^T • (16) 

i=l 

We redefine the ^-coordinate so as to render the metric into the form 

ds 2 = -f(y) 2 dP + dy 2 + a(t, y) 2 j^idx 1 ) 2 . (17) 

i=i 

We call the above bulk-brane configuration the static configuration, since the metric 
components in the (t, ^)-subspace are time-independent, and the branes are not moving. 
Note it was crucial to introduce conformal coordinates before orbifolding. Had we orb- 
ifolded first and moved to conformal coordinates afterwards, the metric would still have 
the form (IT^), but the branes would move through the bulk instead of staying at fixed 
positions. 

What are the constraints for brane matter in the static configuration? If a(t, y) = 
a (t)f(y), the only possible equation of state is that of a cosmo logical constant, p m = —p m - 
Giving up factorizability of a(t, y) does not allow much more freedom. In that case ([13]) 



allows p m to be time- dependent, but because the l.h.s of (14) is time-independent, the 



time-dependent part of 3p m + 2p m must cancel. That corresponds to a time-dependent 
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density of two-dimensional domain walls moving at the speed of light on the brane. So 
the most general equation of state is 



Pm\y) Pm(wall)(t) Pm(vacuum) 



^PmiwalVjit) Pm(vacuum) 



In Section 3 we will discuss constraints for bulk matter. Interestingly, we will find that 
the presence of an ideal fluid in the bulk will force the bulk metric to be that of a static 
configuration. 

Several RS-type solutions containing brane matter with arbitrary equations of state 
have been presented [|], [TT|, [l^, [13], |l^, |2l], [22], ^3], |24| , to mention a few. First of all, it 
is an easier task to find realistic brane cosmologies in a single-br&ne configuration. Then 
one is not addressing the hierarchy problem, and consequently there is more freedom in 
setting up the model. For example, the brane could be moving across the bulk geometry 



or one can take the n(t, y) factor in the metric to be in a non-factorized form. 
For two-brane systems, inclusion of generic matter requires non-static configurations of 
branes. 

We will next briefly discuss the non-static configurations with two branes. 



Non-static configurations. The strong constraints for the brane matter followed be- 
cause we introduced conformal coordinates prior to orbifolding. Let us now see what 
happens if we orbifold in general Gaussian normal coordinates. We start again from the 
general form of the metric @, with c = and K = 0. We then move to the orbifold 
geometry, with parallel branes at the fixed points y\ <y%- Further, we again require that 
6 = and set b = 1. So, we are dealing with the following class of metrics 

3 

ds 2 = -nit, yfdt 2 + dy 2 + a(t, y) 2 ^(rfx^ 2 , (19) 

i=l 

with a compact y-direction y\ < y < y 2 . 

It is known that if n(t, y) takes a non-factorized form, there are single-brane solutions 
allowing generic matter. 

The reference |H| presented a two-brane configuration with generic matter on the TeV 
brane, while satisfying the condition @ up to leading order. However, the fifth dimension 
was required to be static only up to leading order, at higher orders the ^-coefficient in the 
metric may contain time-dependent terms. 



Time-dependence of b. We discussed above that if one assumes b to be fixed, whether 
to avoid inducing time-dependence in particle masses or for some other reason, and have 
found that one easily ends up with overly constrained brane matter. 
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Perhaps the most viable option for including brane matter is to allow for a time- 
dependent b. It has been argued that time dependence in b leads to time- varying particle 
masses [pj. Then any change in b must be negligible by the time of nucleosynthesis so as 
not to contradict observation. As an aside, we note that these arguments have depended 
on a specific form of the metric, and that it is not obvious that a time-dependent b would 
in general lead to such time-dependence of n and a as to affect particle masses. 

We should thus consider the possibility that b is almost stabilized but varies in time 
slowly enough not to conflict with observation, either asymptotically approaching a con- 
stant value or oscillating about a minimum. A b varying slowly should allow the inclusion 
of realistic brane matter. The time-independent b can then be considered an approxima- 
tion, and perhaps an asymptotic limit, of this scenario. Treating time-dependent brane 
matter as a perturbation against a background of bulk matter (including possibly a cos- 
mological constant) and brane cosmological constants is supported by the fact that energy 
density of the universe from the time of nucleosynthesis onwards is quite small in natural 
units, even if the scale of the five-dimensional gravitation is in the TeV range. Also, it 
seems in general more plausible that dynamical evolution would lead to a b with weak 
time-dependence rather than a completely fixed b. 

We now move on to discuss constraints for bulk matter. Note that we will continue 
to take b to be constant, and assume that a weak time-dependence (along with brane 
matter) could be included as a perturbation. This line of thought has been pursued in 



3 The Einstein Equation with Ideal Fluid 



3.1 The stress tensor 

After general investigations of the metric and brane matter, we now proceed to study the 
particular case of a bulk with ideal fluid and a cosmological constant. 

Bulk matter, either in the form of scalar fields or an ideal fluid (which can sometimes be 
interpreted as a scalar field, and vice versa) has been considered in many papers. However, 
the scalar field studies, of which we mention only a few, have concentrated on fixing the 
size of the fifth dimension |4], |^, |(|, inflation |18| |, the cosmological constant problem |20 



singularities and the AdS/CFT-correspondence fI9 |, or on more general aspects of the 
formalism |J, not on obtaining cosmological solutions. In papers of a more cosmological 
nature, often only the yy- component of the bulk stress tensor has been allowed to deviate 
from a cosmological constant []7|, O, [12], |1^, and in any case the yy- component has 
been taken to be different from the afaf-components |TI], [11], |TB[Q, so that the bulk matter 



L6[ is particularly interesting in that the yy- and x^-components are initially assumed to be different 
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cannot be interpreted as an ideal fluid. Furthermore, complete and explicit cosmological 
solutions for the case of a bulk stress tensor with non-trivial tt- and x*x*-components are 
rarely presented; |28j mentions one in passing. We will now present one such solution. 



We take the branes and the bulk to contain some ideal fluids of unspecified nature. 
For the moment, we do not make any assumptions about the metric. It is clearest to 
introduce a local orthonormal frame to find the form of the stress tensor. We introduce 
coordinates (x A ) = (t,x l ,y) such that locally the five- dimensional line element takes the 
form of the 5-d Minkowski metric, 

3 

ds 2 = -dt 2 + ^(d^) 2 + dy 2 . (20) 

8=1 

In the local orthonormal frame, the stress tensor for brane (bulk) ideal fluids must be 
homogeneous and isotropic in the three (four) spatial dimensions on the branes (in the 
bulk). We thus have 

T A 6 = J2 S (y -£m)diag(-p m ,p m ,p m ,p m ,0) 

m=l,2 

+diag(-p - A, p - A, p - A, p - A, p - A) . (21) 

In the above, p m ,p m are the energy densities and pressures of ideal fluid on the two branes 
located at 

ryi rV2 
y x = / dy'b(t,y') ; y 2 = / dy'b(t,y') , (22) 
Jo Jo 

p, p are the energy density and pressure of the bulk ideal fluid and A is the bulk cosmo- 
logical constant which we have for convenience separated out. In particular, note that in 
the local orthonormal frame the pressure of an ideal fluid in the y direction is equal to 
the pressure in the x % directions. We now introduce the assumption of homogeneity and 
isotropy in the directions parallel to the brane. Then the pressures and energy densities 
cannot depend on the coordinates x l , only on the time t and, in the case of the bulk fluid, 
the perpendicular direction y. Thus, expressing the stress tensor in general coordinates 
x A , we have 

T A B = £ 5 y. ^ diag( -p m (t),p m (t) , p m (t) ,p m (t),0) 

m =l,2 \ t - i ym) 

+diag(-p(t, y) - A, p(t, y) - A, p(t, y) - A, p(t, y) - A, p(t, y) — A) . (23) 
The bulk ideal fluid is assumed to satisfy a linear equation of state, 
p(t,y) = wp(t,y). (24) 

but are forced to be equal by the Ansatz used. 
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We could for generality write the bulk ideal fluid as a sum of components with different 
w, but since it turns out that w can take only one value, we prefer not to clutter the 
notation. Also, one could in principle allow the coefficient w to depend on t and y, 
describing a time- and coordinate-dependent (that is, interacting) mixture of ideal fluids. 
However, in this paper we shall assume that w is constant. 

3.2 Bulk ideal fluid implies a static configuration 

We again start from the metric (|3]), with c = and K = 0. Before tackling the Einstein 
equation, we simplify the problem by taking advantage of the conservation law of the 
stress tensor, 

D A T A B = , (25) 
which implies that matter on the branes and in the bulk satisfies the following equations: 



a 

p m + 3(p m + p m ) 

a 



(26) 



y=Vm 



+ + g = (27) 

n! 

p' + (p + p)- = 0. (28) 

Interestingly, ( p8|) implies that it is not possible to have bulk ideal fluid with p = 
(assuming that the warp factor is non-trivial). So, the conservation law of the stress 
tensor already restricts the possible equations of state for the bulk ideal fluid. We will 
shortly see that the full Einstein equation (with the assumption 6 = 0) permits only one 
particular equation of state (aside from the equation of state of a cosmological constant), 
that of the stiff ideal fluid. 

Substituting the equation of state (|24|), we can integrate the equations (|27|) and (|28D 
for the bulk fluid to obtain (for w^O) 

P{t,y) = PoB(t) ™ n[t,y) » 

p(t,y) = p C{y)-^ a{t,y)-^ w n{t,yy^ w \ (29) 

where po is a constant and B and C are some unknown functions, with the powers of w 
introduced for convenience. Excepting the case of a cosmological constant, w = — 1, we 
can rewrite the equations (p9|) as (we introduce the assumption 6 = 0) 

l-\-w . . 1 + w 

P{t,y) = poB{t) - n(t,y) » 

n(t,y) = B(t)C(y)a(t,y) 3w . (30) 



II 



The second equation in (p9|) is a relation between the metric functions n and a. On 
the other hand, the ^-component of the Einstein equation also gives a relation between 
n and a. According to (|), the assumption T ty = implies (with 6 = 0) 



n(t,y) = Ait)' 1 a{t, y) . 
Combining the conservation law relations ( j30|) and the no-flow relation 

l-\-w _ l-\-w 

P\t,y) = PoB{t) « n(t,y) » 



(31) 
we have 



n{t,y) 
a(t,y) 



B(t)C(y)a(t,y) 3w 
A{t)B{t)C{y)a{t,yf w 



(32) 



We can integrate the equation for a in (|3^) and obtain n, a and p in terms of two 
unknown functions of t and two unknown functions of y. The form of these functions can 
be found in the Appendix. 

Let us now consider the Einstein equation (|j) in the bulk, with 6 = (and 6 = 1). 
Since the Einstein equation is local and the metric is assumed to be continuous, the branes 
contribute only to boundary conditions. Hence we can ignore the brane contribution to 



the stress tensor in local bulk calculations. With the relation (plf ), the t^-component of 
the Einstein equation is satisfied trivially. The remaining components of ([D with the 
stress tensor (^) and the equation of state (|24D read, away from the branes: 
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(34) 



Taking linear combinations of the above equations, we obtain the following equivalent 
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3(w + l)p (35) 
{w - l)p - 2A (36) 
. (37) 



12 



Of the above equations, (^) is remarkable in that it does not contain derivatives with 
respect to t. Thus, plugging in n, a and p, which we have expressed in terms of functions 
of t and functions of y, into (|36|), we obtain an algebraic equation involving functions of 
t. The details are again relegated to the Appendix. From this algebraic equation we see 
that n and a have the following factorizable formQ: 

n(t,y) = n (t)f(y) 

a{t,y) = a {t)f{y) . (38) 
Shifting to cosmic time, we arrive at the simple factorizable metric (|9|): 

n(t,y) = f(y) 

a(t,y) = a {t)f(y) . (39) 

So, the presence of bulk ideal fluid leads (given our four assumptions) to the static 
configuration with a simple factorizable metric. 



3.3 Solution of the Einstein equation 

We have shown that when the bulk contains ideal fluid, our four assumptions lead to the 
factorizable metric fl3"9"|). We now complete the solution of the Einstein equation with bulk 
and brane ideal fluids (of course, the brane ideal fluids reduce to cosmological constants, 
as we see from (|13|)). 

With the factorizable metric (|39]), the equations fl35|) to fl3"T|) read 



P = oHzTf ~- U (40) 

6{w + 1) \\a / a ) 

£ = — (4D 

n-w? = -^-|(^) 2 s -£. (42) 

3 ao 3 Veto/ 9 

In the last equation we have introduced a constant C. The l.h.s. of (|41| ) is independent 
of t, while from (fE]) we see that p depends^ on t. Thus, we obtain the result w = lf\. In 
other words, the only allowed equation of state for the bulk ideal fluid is 

p = p. (43) 



7 The special case w = — ^ would yield a non- factorizable metric. However, this requires A = 0, so the 
metric does not have exponential warping and is of no interest as regards the hierarchy problem. 

■ 2 

8 The possibility — = constant would reduce the bulk fluid to a cosmological constant. 
9 This result appeared in [n6| in a different setting, valid only in the case b ^ 0. 
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This result does not depend on the choice K = 0. We postpone discussion of the 
properties of an ideal fluid with this equation of state to section 3, and proceed to find 
the corresponding cosmological solutions. 

With w — 1, the equations (f40|) to fl42|) read 



P 

r 
f 

n - iff 



a 



a 




(44) 

(45) 
(46) 
(47) 



For K 7^ 0, the time dependence of p and the equation for a would change, but the 
equations for the warp factor would remain unchanged. Note that the bulk cosmological 
constant contributes only to the warp factor, not to the cosmological expansion factor a . 
Interestingly, the bulk fluid accumulates at the brane where / has its minimum value, in 



other words at the TeV brane, a result noted in [10] 



Comparison with the FRW model. It is interesting to compare the equations (f44|) 
and ( p|) with the corresponding equations for the K = Friedmann-Robertson- Walker 
model. We take the same matter content for the FRW model: an ideal fluid with the 
equation of state p( 3 ) = p( 3 ) and a cosmological constant A( 3 ) (we use the subscript 3 as 
a reminder that there are only three spatial dimensions). The Einstein equation for the 
FRW model with this matter content can be written as 



+ 2 



P(3) 

a 




3A 



(3) 



(48) 
(49) 



Comparison of ( [44D and ([47]) with (^g) and (|49| ) shows that ao in the RS case is exactly 
the same as a in the FRW case, with C playing the role of a cosmological constant, 
C = 9A( 3 ), and that p equals p( 3 ), save for the factor f~ 2 . These results hold also for 
K ^ 0. So, bulk ideal fluid in the RS model behaves exactly like ideal fluid in the FRW 
model, modulo the warp factor. However, we will shortly see that in the RS case only the 
sign of C has any physical meaning, in contrast to the FRW case, where the magnitude 
of A(3) sets the timescale of the universe. 

We analyze the cases C = and C ^ separately. 
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C = 0. The equations ([44] ) to (|47|) now read 




p = li(-) --)r 2 (50) 

/ 2 

r/-(/0 a = (52) 

/ • \ 2 



(51) 



^ + 2 p) =0. (53) 

CLq \CLq, 



The above equations have the following unique^ solution 



a (t) = (- (54) 



T 



f(y) = f e ±y ^ 2 (55) 

P(t,y) = -t^ftvr 2 - (56) 

We have normalized by ao(0) = 0. The bulk cosmological constant has to be negative, 
A < 0. 

The solution contains two free parameters, r and Jq. A choice of the time parameter 
t corresponds to choosing the unit of time, so that r has no physical meaning, as in the 
FRW case. A choice of the warp parameter f corresponds to choosing the origin of the 
y-coordinate. The metric and the energy density are invariant under the scaling t — > At, 
x % — > Xx\ f — > A _1 / , so that fo, or the placement of the origin, has no physical meaning. 

The C = solution contains no physical degrees of freedom other than the value of 
the bulk cosmological constant: the metric and the bulk energy density (and, as we will 
see, the brane cosmological constants) are fixed once the bulk cosmological constant is 
specified. This suggests that the solution is unstable, possibly collapsing to one of the 
C / solutions when perturbed. 

C In this case the equations ([0]) to ( |4"TD read 

= M r2 (57) 

6 \ V a / a J 
f" A 

7 = "I (58) 

f"S-U'f = (69) 

/ • \ 2 
fl , o / a 



+ 2 - = i- (60) 



10 



Apart from the trivial solution do = 0, which would lead to p = and the original RS model. 
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The above equations have the solutionF] 



, sinV^JlCi t) C<0 f . 

a ° {t) = smh^t) c>0 (61) 

C 



fiv) = /ie^ + -y- e ^ (62) 



fk,v) 3k I sinh- 2 (VCt) C>0 V ; 

We have again used the normalization ao(0) = 0. The equations for / would allow a 
positive (zero) bulk cosmo logical constant, resulting in trigonometric functions (a linear 
function), but in order to solve the hierarchy problem we take A < 0. 

As in the previous case, the parameter a\ has no physical meaning, and can be set to 
unity, while fi corresponds to the choice of the origin of y. Rewriting fi in terms of a 
coordinate value yo, we have 

= I Sin !i 3 iy^:; (64) 

M = K =^{ —vv^rv vo)) c<o 

* " ' ' / |A|/2(y- W) )) C>0 




pM) = ™ vvn.. -y ))^-y\c\ t) c<o 

• °(V|A|/2 (s/-jA))) sinh-VC*) C>0 

We see that the metric and the energy density are invariant under the scaling t — *■ At, 
x l — ► Ax*, C — > A _2 C. In other words, only the sign of C has any physical meaning, the 
magnitude is irrelevant. In particular, C does not introduce a new mass scale into the 
model. 

Unlike in the case C = (and the original RS proposal), the model is not invariant 
under translations of the y-coordinate, so that the choice of origin of y does have a physical 
meaning. The translational invariance is broken by the ^-dependence in the bulk matter: 
y is not a free parameter but is set by the bulk energy density, as we see from ( |66"D . There 
are two physical degrees of freedom: the value of the bulk cosmological constant and the 
bulk energy density. (In addition, there are of course the brane cosmological constants, 
which, as we will see shortly, are also free parameters.) 

In the case C > the function / has a zero at y = yo. The bulk energy density and 
the scalar curvature diverge when / = 0, so the singularity at y — yo is physical, not an 

11 The equation for ao allows de Sitter and anti-de Sitter solutions when C > 0, but these would again 
give p = 0. These solutions were first presented in m% E9|. 
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artifact of the coordinate system we have chosen. In what follows, we simply avoid the 
singularity by constraining the value yo not to lie between the brane positions in the case 
C > 0. 



3.4 Branes and fine-tuning 

Having solved the Einstein equation in the bulk, we now turn to the branes, which, as 
noted earlier, provide boundary conditions. 



Brane cosmological constants. The factorizable metric implies that the branes con- 
tain only cosmological constants. We switch to the notation p m = —p m =: A m . Inserting 
the brane contribution to the stress tensor from fl23| ) into the Einstein equation ([|), we 
get 



/ 



(67) 



y=yv 



where yi,y2 (y2 > Hi) are the orbifold fixed points, where the branes are located^. As 
before, the notation [/'] refers to the discontinuity of /'. Writing the discontinuity in 
terms of the continuous part f , we have 



f 

9_ 
/ 
f 

r\ J C 
f 



y=yi 



y=V2 



-Ai 



-A, . 



(68) 



For C = (and the original RS proposal) we have f(y) 
gives the result Aj : 
have instead 



foe' 



Then (Eg 



-A, 



2|A|, an unfortunate fine-tuning. However, for C < we 



Ax = - v/2^tanh( yJ\£\/2 {y t - y )) 
A 2 = ^jA|taDh(^/jAj72(j/ 2 -j/o)) ; 



(69) 



for C > the hyperbolic tangent is replaced by a hyperbolic cotangent^. The above 
equations for the case C < were found in in a slightly different setting. The values 
of hyperbolic tangent lie in the interval from -1 to +1 (excluding ±1), and the values of 
hyperbolic cotangent stretch from minus infinity to plus infinity, excluding the interval 

12 In the case C > we demand tjq < y\ or yo > y% to avoid having a singularity. 
13 Note that either set of equations fixes y% — y%, the size of the fifth dimension. 
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from -1 to +1. Thus, the ratio of the brane cosmo logical constants to the square- root of 
the bulk cosmological constant determines the value of C: 



C < for -^y=<l 
C = for -^y= = l 

C > for -44= > 1 • (70) 
2|A| 

In the present model, the Einstein equation implies no fine-tuning of parameters, unlike 
in the original RS proposal. For a given value of the bulk cosmological constant A, the 
absolute values of both brane cosmological constants have to be smaller than, bigger than 



or equal to y2|A|, but are otherwise unrestricted by the Einstein equation. With two 
equations to satisfy and two constants, y\ — y and y 2 —yo, at our disposal, no fine-tuning 
is needed to obtain a solution to (|68|). 

The signs of the brane cosmological constants are opposite for C = and C > 0. 
This can also be the case for C < 0. However, for C < it is also possible for both 
brane cosmological constants to be positive: this requires the value of yo to lie between 
the brane positions, y\ < y < y 2 . Furthermore, we have the interesting possibility of 
a zero cosmological constant on one brane and a positive cosmological constant on the 
other: this requires one of the branes to be placed at the point where / has its minimum, 
Ho — Ui or Ho — Vi- (It is n °t possible to have two branes with non-positive cosmological 
constants.) Thus, there is no need for negative energy densities on the branes, essentially 
because the warp factor is not a monotonical function, but has a minimum. This was first 



noticed in |T_7[ and further emphasized in 0. Of course, the bulk cosmological constant 



still has to be negative to obtain an exponential warp factor. 

The hierarchy problem. The solution to the hierarchy problem introduces a new 
equation, namely 

f(yi) . 10 ie (71) 



For C < 0, the above equation reads 

cosh(y^72 (m - y )) 
cosh( J\X\/2 (y 2 - y )) 



10 16 ; (72) 
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for C > the hyperbolic cosine is replaced by a hyperbolic sine. Combining ([72]) with 
(|69|), we can write the hierarchy condition ([n]) as the following equation, valid for all 
values of C: 



A? „ / A 



KT* - 1 . (73) 

2|A| V2|A| J 

We see that the only way to avoid introducing unnaturally large numbers is to fine-tune 
the cosmological constants, A^ = A% = 2|A|, corresponding to C = 0. 

The above result is not dependent on the presence of bulk ideal fluid: if we put p = 0, 
our solutions disappear, but we get new solutions with the same problems. For C = we 
get a constant ao (the original RS solution) and for C > we get de Sitter and anti-de 



Sitter solutions, found in [[H], |29j; the case C < becomes disallowed. The warp factor is 
unchanged by the absence of bulk ideal fluid (except that y becomes a free parameter), 
so the condition ( |73|) is also unchanged. This means that the fine-tuning problem is not 
inherent to the Einstein equation. Specifically, it is not due to the condition b = 0, as 
sometimes claimed It is rather a rephrasing of the hierarchy problem. In this sense 
the RS model does not provide quite a satisfactory solution to the hierarchy problem. 



4 Discussion 



As is well known, the RS solution is quite precarious; bulk and brane sources for gravity 
must be chosen very carefully for the RS solution to emerge from the Einstein equation. 
However, it turns out that in two-brane systems one easily end up with severe constraints 
for the sources; this fact is not apparent from the literature. We hope that our results 
help to identify the potential problem sources in trying introduce realistic matter into 
the bulk or branes, while preserving the main features of the RS scenario: the hierarchy 
generation and the localization of gravity. The need to include generic brane matter in 
particular creates problems in trying to interpret our universe as a brane world. 

In Section 2 we found that under minimal assumptions on the RS model, a stabilized 
radion and zero thickness branes, one cannot introduce realistic brane matter^] in static 
two-brane configurations. The only allowed equation of state, other than that of a cos- 
mological constant, corresponds to a fluid of two dimensional domain walls moving at the 
speed of light. Even that requires giving up factorizability of the bulk metric. Further re- 
laxation of the constraints for brane matter requires a time-dependent and non-factorized 
warp factor, and most likely a slow stabilization of the radion as the branes settle to the 
final fixed proper distance. This idea is supported by the result of Section 3, which shows 

1 In this context, we mean homogeneous and isotropic matter with equations of state corresponding 
to a combination of dust and radiation. 
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that when the bulk is filled with ideal fluid and the radion is completely stabilized, the 
warp factor is necessarily factorizable. The constraint for the equation of state of the bulk 
ideal fluid also points in this direction. 

The dynamical stabilization of the radion degree of freedom can be expected to be 
highly nontrivial. Nevertheless, unless for some reason the universe started out with 
all the degrees of freedom in the ground state, some dynamics should be expected. In 
the context of cosmology, a natural initial condition for the radion could be any value 
compatible with the uncertainty principle, as in chaotic inflation. The radion potential 
energy must then be released in some way. It is conceivable that the extra radion energy 
is dissipated into purely gravitational degrees of freedom, but there is no guarantee that a 
RS-type solution would result; this remains to be studied. Dissipation into fields living on 
the brane has been studied in || . The third possibility, radion decaying into bulk degrees 
of freedom, assumed in this paper, is possible but highly constrained, as we have shown. 
The only admissible equation of state for ideal fluid in the bulk was found to be p = p , 
representing the so-called stiff ideal fluid. Here " stiff" reflects the fact that the velocity of 
sound in the fluid is equal to the velocity of light. Concretely, such a fluid corresponds to 
a classical free massless, coherent scalar field (not to be confused with massless radiation), 
something that was considered in a static setting in . Further, the massless bulk scalar 
field should not be confused with the massive Goldberger-Wise bulk scalar field. Whether 
such a field could be coupled to the radion is beyond the scope of the present study. 

The no-go flavor of our results is meant to clarify what are the most promising direc- 
tions for modifying the two-brane models so that other equations of state for the brane or 
bulk ideal fluids can be allowed. In the case of the bulk, a natural guess would be to start 
with a different Ansatz for the bulk geometry. This corresponds to a different assumption 
about the spacetime symmetries. For example, if the bulk geometry would correspond 
to an AdS black hole, one might expect the equation of state for massless radiation to 
become allowed. 
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Appendix: Bulk Ideal Fluid and Factorization 



We present here details of the calculation of section 3.2 which shows that the presence of 
bulk ideal fluid leads to a factorizable metric. The starting point was the equations (R2T) : 



P{t,y) = PoB{t) - n(t,y) » 
n(t,y) = B(t)C(y)a(t,yf w 
a(t,y) = A{t)B{t)C{y)a{t,yf w 



(74) 



Integrating the equation for a, we obtain the following expressions for n, a and p: 



1 

w = — 
3 



w 7^ 



p(t,y) 

n{t,y) 
a(t,y) 

P(t,y) 
n(t,y) 
a(t,y) 



p C{yY A e -^C{y)E(t)+D{y)) 

B{t) C(y) e c ^ E ^ +D ^ 

e C(y)E(t)+D(y) ? 

Po C(y)- 1 ^(C(y)E(t) + D(y))-^ 
B{t)C{y) (C(y)E(t) + D(y))^ 
(C{y)E(t)+D{y))T^ . 



l + w 



(75) 



In the above expression, D(y) is a new unknown function, and E(t) is proportional to 
the integral of A(t)B(t). We have now separated the t- and y-dependence of the metric. 
The linear combination of components of the Einstein equation which does not involve 
derivatives with respect to t is (BB): 



a" n" 
3— + — 

a n 



[w 



l)p-2A . 



Now we simply insert (|75|) into (76). Let us consider the cases w 
separately. 



(76) 

| and w 7^ | 



w = |. In this case ( |76D reduces to 



/2 s 



AC ,2 E(t) 2 + Ud" + SC'D' + 2^- j E(t) 



All fil 

AD" + AD' 2 + — + 2=^D' 



A 



C 



w 



l)p C' A e-^ CE ^-2A. 



It is impossible for the above equation to be satisfied at all times t unless E 
C = 0, both of which correspond to a = 0, and thus static cosmology. 



(77) 



or 
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w 7^ ~. Now (|76|) can be written as 



Gw(l + 3w) (D> C'V 
3(1 + ^) /D" C" 6w C ( D' C'\ , 



+ 




3w C 



/2 X 



1 - 3w \ C 1 - 3w C 2 J 
= (w - 1) poC-'^r (CE(t) + D)- 3 ^t -2A. (78) 

The above equation is less transparent than ([771) . Unless u> = — |, there are three 
powers of CE + D on the l.h.s., so by equating the coefficients of the powers of CE + D we 
obtain three equations. These equations can only be satisfied if D = (or D oc C, which 
amounts to the same thing by a redefinition of C). From (|75|) we see that this corresponds 
to a factorizable metric. This solution also requires w = 1, though for reasons of clarity 
we preferred to present this result at a later stage in section 3. 

As for the special case w = — ~, the above equation then does have solutions with a 
non- factorizable metric. However, a substitution to the full Einstein equation (|35|)-(|37l) 
shows that this requires A = 0, and as a result the metric does not have an exponential 
warp factor and is of no interest as far as the hierarchy problem is concerned. 
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